Field Theory of Critical Behaviour in Driven 
Diffusive Systems with Quenched Disorder 

V. Becker and H.K. Janssen 
Institut fiir Theoretische Physik III 
Heinrich-Heine-Universitat Diisseldorf 
D-40225 Diisseldorf, Germany 

March 22, 1998 



We present a field theoretic renormalization group study for the criti- 
cal behaviour of a uniformly driven diffusive system with quenched disorder, 
which is modelled by different kinds of potential barriers between sites. Due 
to their symmetry properties, these different realizations of the random po- 
tential barriers lead to three different models for the phase transition to 
transverse order and to one model for the phase transition to longitudinal 
order all belonging to distinct universality classes. In these four models that 
have different upper critical dimensions dc we find the critical scaling be- 
haviour of the vertex functions in spatial dimensions d < dc. Its deviation 
from purely diffusive behaviour is characterized by the anomaly-exponent r] 
that we calculate at first and second order, respectively m. e = dc — d. In 
each model rj turns out to be positive which means superdiffusive spread of 
density fluctuations in the driving force direction. 
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1 Introduction 



For more than a decade the long-time and critical behaviour of diffusive sys- 
tems subjected to a driving force has attracted considerable interest. This is 
mainly caused by the richness of their highly nontrivial features which gener- 
ically result from the fact that, due to the driving force, in general even the 
steady states of these systems are far from thermal equilibrium, especially 
in the case of very strong driving forces. Further, driven diffusive systems 
might be suitable models for fast ionic conductors which was first suggested 
by Katz, Lebowitz, and Spohn A review of the different investigations 

on driven diffusive systems and their relations to other non-equilibrium sys- 
tems is given by Schmittmann and Zia 0. 

We are interested in the diffusive motion of uniformly driven particles with 
short range attractive interactions and hardcore repulsion. With such an 
interaction driven diffusive systems show two kinds of phase transitions. At 
the transverse (with respect to the driving force E) phase transition the sys- 
tem changes from a disordered state to an ordered one where the system 
is ordered in the transverse direction and remains disordered in the longi- 
tudinal (parallel to the driving force) direction. The typical configurations 
are strips of high- and low-density phase, arranged parallel to the driving 
force E (Fig. |l](a)). At the longitudinal phase transition the systems changes 
from a disordered state to a state where the system is ordered in the longi- 
tudinal direction and remains disordered in the transverse directions. Here 
the typical configurations are domains ("pancakes") of different phase with 
interfaces perpendicular to the driving force (Fig. |l|(b)). In preceding pa- 
pers the long-time and critical behaviour of driven diffusion in an ordered 
medium and the long-time behaviour of driven diffusion in a medium 
with quenched disorder [0 have been investigated by renormalized field the- 
ory. In the present paper we study the effect of quenched disorder on the two 
kinds of phase transitions in driven diffusive systems. Quenched disorder is 
important in real systems, as it models impurities and defects. We stress that 
this work completes the investigations of a whole model class that is graph- 
ically shown in the summary (Fig. |^). Note that the effects of quenched 
disorder in randomly driven diffusive systems were recently analyzed 
There the driving force is itself a locally random variable with respect to its 
amplitude and sign, whereas here the driving force is uniform. 
The paper is organized as follows: In Section II we set up a Langevin descrip- 
tion of our system, directly on a mesoscopic length scale using conservation 
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laws and symmetry arguments. The quenched disorder is modelled by ran- 
dom potential barriers between sites, and the symmetry properties of the 
random potential play a crucial role. Different realizations of the random 
potential are possible and are argued to lead to different kinds of noise. In 
Section III we study the transverse phase transition in the convenient for- 
malism of the dynamic functional. In Section IV the longitudinal phase 
transition is analyzed. Section V contains the main results and a graphical 
overview of the entire model class. 
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2 Model Building 



The configurations of a driven diffusive system witli liomogenous driving force 
are fully characterized by a single conserved scalar field, i.e., the local particle 
density n{r,t). As the particle number is conserved the order parameter for 
both phase transitions is the deviation of the actual density n(r, t) from its 
uniform average no: 

s{r,t) = n{r,t) - uq . (1) 
This fluctuating variable satisfies a continuity equation 

s + V-j = (2) 

where the current density j consists of a deterministic and a random part. 
The resulting stochastic differential equation is a Langevin equation. First we 
model the deterministic part. Due to the anisotropy of the system caused by 
the driving force, j shows different symmetries longitudinal and transverse 
to its direction. Henceforth the indices || and ± denote the spatial direc- 
tion longitudinal to the driving force and the {d — 1) -dimensional subspace 
transverse to it, respectively. The transverse current is according to 

j± = -AV_L//_L (3) 

caused by a chemical potential with A being a kinetic coefficient. As 
no direction is selected in the {d — 1) -dimensional subspace the transverse 
current jj. is isotropic and therefore a vector of the form 

j± = Vx/(s,A^,V||) (4) 

where the Laplacean Aj^ as argument denotes an even number of V_L-operators 
in every term of /. An expansion of / with respect to s and gradient- 
operators yields, up to higher order terms 

j± = V±(s + s2 + Axs + V||S + A||s) , (5) 

where here and in the following equations coefficients are suppressed for the 
sake of simplicity. 

The longitudinal current, however, possesses no symmetry, due to the driving 
force. As a scalar it can be written as 

i|l=^(s,V||,A^) (6) 
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where the isotropy of the transverse subspace is again taken into account by 
the even number of V_L-operators. Its expansion up to higher order terms 
reads 

j|l = c + s + + V||S + Vjs + Vfs + V||S^ + A±s + V\\A±s . (7) 

While some terms in this equation originate in a chemical potential /xy, ac- 
cording to 



others are due to the driving force E, which at least produces the terms 
proportional to c, s, s^, and A±s and which, in principle, also contributes 
to all terms of the longitudinal current in (J^). Here, the constant c is the 
homogenous part of the current, and is the first nonlinear term of the 
longitudinal current and so the leading nonlinearity of the problem. 
Note that a repeated coarse graining of a microscopic model typically gener- 
ates anisotropic transport coefficients. Building here a model directly on a 
mesoscopic length scale we must take that into consideration. Thus, although 
containing the same terms, /iy and fi± in general have different coefficients, 
due to the anisotropy. From a technical point of view anisotropic transport 
coefficients are required to make the model renormalizable. 
Following microscopic models and simulations [[ll, p|, p!0|-|T2| we restrict our- 
selves to such models which additionally hold a CP (charge and parity)- 
symmetry, i.e. in which the Langevin equation is invariant under the trans- 
formation 



In microscopic lattice models which also form the basis for Monte Carlo sim- 
ulations of driven diffusive systems the PC-symmetry corresponds to the 
particle-hole symmetry in the half-filled system. 

In models with PC-symmetry all terms from and vanish whose sum 
of s- and Vy-factors in the Langevin equation is even. 
The stochastic part of the Langevin equation ( = —V ■ Jl reflects a random 
current that summarizes the fast microscopic degrees of freedom (local in 
time in our Markovian continuum description) and the effects of the quenched 
disorder of the medium. After writing down the general form of the Langevin 
equation we show how to model the various possibilities of quenched disorder 
and which parts of the noise are relevant in the renormalization group sense. 
Note that the noise considered here is conserving, i.e., it satisfles the continu- 
ity equation. We mention that we have also analyzed driven diffusive systems 




(8) 




(9) 
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with nonconserving noise that is caused by random particle sources |TB, H . 
Taking anisotropic transport coefficients into account we obtain the Langevin 
equation 



s = A[A^(rx-/txA^)+pA||(r||-K,||A||)-/tAxA||]s + -A^V||S^-V-jL • (10) 

This constitutes the fundamental equation for driven diffusive systems with 
homogenous driving force, attractive interactions, CP-symmetry and con- 
serving noise. It contains all terms which are relevant in the renormaliza- 
tion group sense for the two phase transitions and the noncritical disordered 
phase, but in each case it still contains irrelevant terms which must be elim- 
inated by a dimensional analysis. 

We proceed by investigating the influence of the quenched disorder in detail. 
In a microscopic driven lattice gas model the quenched disorder is modelled 
by random potential barriers between the sites. There are three possible re- 
alizations of their randomness which is depicted for one dimensional systems 
in Fig. ^ 

I The particles are in randomly deep potential valleys that are separated 
by randomly high potential mountains (Fig. ^ (a)). 

II The potential valleys are randomly deep, the potential mountains are 
equally high (Fig. |^ (b)). 

Ill The potential mountains are randomly high, but the potential valleys 
are equally deep (Fig. || (c)). 

The homogeneous driving force tilts this landscape in the longitudinal direc- 
tion by an angle that depends on its strength. Thus in this direction the 
symmetry of the random potential in the realizations II and III is broken, 
whereas the symmetries in the transverse subspace are not affected. Hence 
in the driving force direction the random potential is unsymmetric in all re- 
alizations. 

In the microscopic lattice gas model particles can only jump to neighbouring 
unoccupied sites with jump rates that depend on the external driving force, 
on the energetic situation of the particles due to their attractive interaction, 
and on the locally random height of the potential barriers between the sites. 
Performing a continuum limit of the microscopic model one can show that, 
first, the main effect of the quenched disorder results in a time-independent 
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random current with zero mean and Gaussian fluctuations and that, second, 
the three different reahzations of the random potential produce three differ- 
ent types of such a random current Cd(r)- 

In the unsymmetric case (realization I) the correlations of the random current 
are given by 

<Cd(r)> = (11) 
< Crf(r) ® Cd(r') > = 2A25(r-r')[(Te|| ®e|| +7(1 -e|| ®e||)] . 

Thus the correlations of the noise force V ■ (d read 

<V-Cd(r)> = (12) 
< V ■ 0(r)V ■ C(r') > = -2A2(7Ax + (TA||)5(r-r') . 

By a suitable scale change of s the kinetic coefficient A in equation(p!2[) is the 



same as in equation(|lOD. The anisotropy of the system due to the driving 

force is taken into account here by the transport coefficients 7 and a. 

In realization II a particle at a given site sees equally high potential barriers 

in all transverse directions. As a consequence of this symmetry the random 

current in the transverse subspace vanishes in the lowest order and starts 

with a local gradient. The correlations of the random current are here given 

by 

< V ■ C(r)V ■ C(r') > = -2A2(-aAi + fTA||)(5(r - r') . (13) 

As the symmetry of the random potential is broken by the driving force in 
the longitudinal direction the longitudinal random current behaves as in the 
unsymmetric case. 

In realization III a potential barrier looks the same from both sides in the 
transverse subspace. This is why the transverse random current vanishes 
totally. Being only longitudinal the random current has the correlations 

< V ■ C(r)V ■ C(r') > = -2AVA||5(r - r') . (14) 

Since the continuum limit of the microscopic model is not rigorous, it might 
be possible that even in the symmetric realizations II and III transverse noise 
terms with a ^ and 7 7^ 0, respectively, are generated by the coarse grain- 
ing procedure. This is due to the fact that the symmetry of the random 
potential is microscopic and concerns only a single site or bond. As we will 
see later each type of noise nevertheless leads to different critical behaviour 
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governed by different stable fixed points with finite regions of attraction. We 
defer tlie discussion of tlie possible additional transverse noise terms until 
the end of the Sections 3.2 and 3.3, respectively. 

By dimensional considerations (power counting) one can prove that each of 
the equations (|12|) - ([T^) contains all relevant noise terms in the renormal- 
ization group sense for the respective type of the random potential. Also 
deviations from the local Gaussian nature of the random current ( are irrel- 
evant for the critical and long-time properties. Even the microscopic part 
(m of the random current summarizing the fast microscopic degrees of free- 
dom is irrelevant if (d 7^ 0. It should be remarked that (m plays the role of 
a dangerous irrelevant field for the calculation of correlation functions with 
frequencies u ^ 0. Without Cm correlation functions show only a "central 
peak" at u; = 0. 



8 



3 Transverse Phase Transition 



In the following all models are analyzed with the help of renormalized field 
theory. This method has successfully been applied to all driven diffusive 
models investigated so far. 

The three different realizations of the random potential described in the 
last section actually lead to three different models for the transverse phase 
transition and have to be treated separately. 



3.1 Unsymmetric Random Potential 



The model for a driven diffusive system with frozen random unsymmetric 
potential is based on the equations (|TUp and ([T^). 

To set up a renormalized field theory, it is convenient to recast the model in 
terms of a dynamic functional [jT^|-|2^ 

J[s,s\=j(fr^jdt[s (s + A(A^(k_lA_l - r^) +pA||(K||A|| - ry) + k^^/\\\] 



„2l 
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A JdtV^s 


— CT 


A jdtV\\s 





(15) 



where s(r,t) is a Martin-Siggia-Rose |^T|] response field. Correlation and 
response functions can now be expressed as functional averages with weight 
exp(-J). 

But for the description of the transverse phase transition this dynamic func- 
tional still contains irrelevant terms in the renormalization group sense. Now 
these terms are determined by a dimensional analysis. 

The transverse phase transition is characterized by finite ry and t±_ —>■ 0. We 
introduce a scale /x^ for small t±. Then fi~^ is a convenient length scale. 
Since t± tends to at the transverse phase transition, the leading term in 
the transverse direction is proportional to sXAj_k,j_Aj_s, whereas in the lon- 
gitudinal direction the leading gradient term is proportional to sXpA\\T\\s. 
The comparison of the leading gradient terms demonstrates that Ay scales 
as A^. For longitudinal and transverse length scales, this implies 

r± ~ yU~^ ry ~ yU^^ . (16) 

Since the dynamic functional is dimensionless the dimensions of fields and 
coupling constants are 
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d+7 



9 



9-ci 



7 ~ /i 5( ~ /i 2 , (17) 

i.e. the coupling constants Ky, k, and a are irrelevant in the renormalization 
group sense. The irrelevancy of a indicates that for the transverse phase 
transition in a unsymmetric random potential the noise in the driving force 
direction is irrelevant. From the dimension of the nonlinear coupling constant 
g we recognize 

4 = 9 (18) 

as upper critical dimension of this model, above which g is irrelevant und 
below which g is relevant. The finite ry and the transverse coupling constants 
are absorbed into lengths and fields by a suitable scale change. Thus, the 
appropriate dynamic functional to describe the transverse phase transition 
in a driven diffusive system with an unsymmetric random potential is given 
by 



~ss + A5(Ax(Ax - r^) - pAj|)s + 1a^(V||5)s 

A IdtV I s 




The dynamic functional has the following symmetries. The isotropy in the 
transverse subspace and the CP-symmetry that reads here 

rj| — > — ry s —s s ^ —s (20) 

were directly integrated in the model building. After averaging over the 
quenched disorder, J' exhibits translational symmetry in space and time. 
The invariance of under the longitudinal scale transformation 

ry — >■ (3r\\ r± — >■ r± (21) 

s — > f3~2s s — >• (3~2s 
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is of great importance, because the parameter combination g'^p~'i is found as 
appropriate variable of the model being invariant under this transformation. 
In contrast to the corresponding model without quenched disorder the 
dynamic functional is here not invariant under a Galilei transformation. 
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To study the critical properties of the transverse phase transition we apply 



standard renormalization group methods |jT9|,|2y, We use dimensional 
regularization in d = 9 — e followed by minimal subtraction. The one-line- 
irreducible vertex functions with n S-legs and n s-legs at wavevectors {q} 
and frequencies {uj} will be denoted by r^j^^dq}, {oj}). Taking into consid- 
eration the causality properties of the theory we only find Fi i and Fi 2 prim- 
itively divergent. The nontrivial diagrams contributing to Fj^ ,i({q}, {a;}) 
carry at least a factor q^, because the interaction vertex is conserving. The 
primitive divergences are multiplicatively absorbed in a redefinition of the 
parameters 

p p = Zpp , g g = H^ZuU . (22) 

Here, in contrast to the ordered problem the coupling constant g has 
to be renormalized, due to the loss of Galilean invariance in the disordered 
problem. 

The elements of our perturbation expansion are the Gaussian propagator 

iu; + A[q]_(q]_ + rj_) + pgyj 



the Gaussian correlator 



u;2 + A2[q^(q^ + r^)+pg2 



and the conserving vertex f(q) = —i\gq\\. Their graphical representation is 
shown in the Appendix in fig.(^). 

As the correlator of the theory Cq^^^ is proportional to 5{(jj) (i.e. Cq(t) is 
independent of time) and any loop of a diagram contains at least one corre- 
lator because of the structure of the vertex with two s-legs and one s-leg and 
because of causality the integration over internal frequencies cannot generate 
any divergences. Thus the divergent parts of the vertex functions are only in 
the {uo = 0)-parts. To facilitate their calculation it is convenient to split the 
order parameter 

s(r,t)=^(r)+s'(r,t) (25) 

into a time independent (p{r) and a time dependent s'{r,t) that contains 
no (a; = 0)-parts. Since the noise is independent of time s'{r,t) relaxes 
deterministically. In the limit of long times we obtain 

s(r,t ^ oo) = v?(r) . (26) 
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Together with 

<f{r) = A y dts{r, t) (27) 
the dynamic functional reduces to a quasi-static (frozen) Hamiltonian 

n[^,<f]^ Jd''r^^<f{A4A^-Tl_)-pA\^)if+^g{V\\<f)ip^ + <fA^<f^ . (28) 

The frozen Hamiltonian generates all the zero-frequency parts of the vertex 
functions if causality is included in the graphical rules of perturbation theory. 
In the quasi-static model the Gaussian propagator and correlator are only 
dependent on wavevectors and read 

^ qi(qi + rJ+p?J ^^^^ 



[qi(qi + rx)+pgf]^ 

In the quasi-static model we calculate the primitively divergent vertex func- 
tions in two-loop order. While the one-loop calculation is easy to perform 
analytically, technical difficulties have to be overcome for the two-loop di- 
agrams with mixed q^-propagators and -correlators. These two-loop dia- 
grams have here been solved by a new technique where a special inverse 
Mellin transformation is used to factorize the denominators of some propaga- 
tors and correlators, respectively. Thus, the e^-poles and the simple e-poles 
can be extracted. Their coefficients are then given by parameter integrations 
over paths in the complex plane. Whereas the coefficients of the e^-poles can 
be analytically calculated, those of the simple e-poles have to be computed 
numerically. The details of the two-loop calculation are shown in the Ap- 
pendix. 

In a two-loop calculation using dimensional regularization we obtain the 
vertex functions in an expansion in e = dc — d and q 

f M(q) =qi(qi + r^) + PQ^ + ' ^^^^^^(l - 0-2158e) 

f i,2(q) - ^l — iq\\r^ + " 0-2154e) , (30) 

D e p2 o 6 p 

where bare unrenormalized parameters are indicated by a superscript " ° " 
above the symbol. — j^jgOd-ir{l + |)r(^)r(|) is a suitably chosen 
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constant with Od-i being the surface of the {d — 1) -dimensional unit sphere 
and T{z) denoting Euler's F-function. The Z-factors defined in equation (p2D 
are in minimal subtraction 

Z, = l--AA-f^^?^(l + 0.2158e) + O(«6) (31) 
3e pi de-' p-^ 

Zu = 1+ A,^ + —A'^^il + 0.215Ae) + O{u'). 
be p2 8e^ p-^ 

We recognize that the perturbation expansion is organized in powers of the 
dimensionless renormalized parameter combination v := A^u'^p~2 that was 
already found as invariant variable under the longitudinal scale transforma- 
tion (equation(pID). 

With the renormalizations at hand we are in a position to determine the 
critical behaviour of the vertex functions. We use the fact that the unrenor- 
malized theory is independent of the momentum scale p. This leads to the 
renormalization group equation 

rn,n({gii,q±},T^±,^^,p,Ai) = . (32) 

The Wilson parameter functions, being only dependent on v, are given by 

v{e -^v- 0.252?? + 0{v^)) (33) 

--V - 0.0959^2 ^ Qf^^s^ ^ 
3 

where the derivatives are calculated at fixed bare parameters. As a linear 
partial differential equation the renormalization group equation is solvable 
by the method of characteristics with the result 

rn,n({g||,q±},rj_,t;,p,/i) = rfi,„({g||,q_L},r^, ?)(/), p(/),/i(0) • (34) 
The trajectories v{l), p{l), and p{l) are solutions of the flow equations 

ij^vii) = pMi)) ij^ inp(/) = am ijf^i) = m m 

with the flow parameter I and the initial conditions 

v{l = l)=v p{l = l) = p p(l = l) = p. (36) 



P ^ p(^^ p^ 
" dv dp dp 



dv 
^d-p 
dlnp 



dp 
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The third flow equation obviously has the solution = fil. 
With the help of characteristics the critical asymptotic region of small t± 
and q can be mapped onto uncritical regions. In the scaling limit / ^ 1 
corresponding to 



< 1 



< 1 



< 1 



(37) 



the flow of v{l) is controlled by stable zeroes of Py. For e > 0, there is a 
nontrivial infrared stable fixed point 
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(38) 



The other fixed point f * = being Gaussian is stable only if e < 0. 

From the nontrivial fixed point value C('^^*) we define the anomaly exponent 

r] := — |C(^*) that is to two-loop order 



(39) 



V = ^e{l- 0.302 - + 0{e')). 



At the fixed point the solution of the second fiow equation of equation 
is 

m = pi-^' ■ (40) 

The results found for the vertex functions in the quasi-static model can be 
directly transferred into the dynamic model because the time scale A needs 
no renormalization. 

The longitudinal length scale transformation according to equation (PT|), di- 
mensional analysis and the renormalization group equation are now combined 
to derive the asymptotic critical scaling form of the vertex functions in the 
dynamic model at the transverse phase transition 

r«,n({g||,q±,cj},T-j_,t;^, A,p,/i) = 

This equation implies that, below the upper critical dimension rfc = 9, anoma- 
lous scaling behaviour only occurs in the direction of the driving force and is 
completely characterized by the anomaly exponent ?7 from ( p9| ) that is posi- 
tive. These results are in analogy to what was found in the driven diffusive 
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systems investigated so far [P,|5|,[7|,P!B[]. 

To illustrate the importance of the positive anomaly exponent rj we especially 
investigate the density response function x(q, t) which is the Fourier trans- 
form of T{\(q^,u). Choosing the flow parameter / = cu* <^ 1, we derive the 
scaling form of the density response function x{^.jt) •=< Q, 0) > 

from equation (p^: 

X(q,t) = . (42) 

From that we conclude that for long times typical longitudinal length-squares 
scale with time as 

< rJ > ~ . (43) 

The positivity oi r] = ^e(l — 0.302 g + 0(e^)) means that in systems with 
spatial dimensions d < = 9 fluctuations spread faster than diffusively 
in the driving force direction. This superdiffusion was also observed in all 
driven diffusive systems analyzed up to now P,|5|,[f|,p!3[]. 
A comparison with the corresponding model for the transverse phase transi- 
tion without quenched disorder, having an upper critical dimension dc = 5, 
shows that the phenomenon of superdiffusion occurs here over a much greater 
dimensional range. As the system without quenched disorder behaves nor- 
mally diffusive for c/ > 5 it follows that quenched disorder is the reason for 
superdiffusion in the dimensional interval 5 < < 9. 

The analogous comparison of the models with and without disorder for the 
noncritical region came to an analogous result. 

Another consequence from the scaling form of the density response function 
is that typical transverse length-squares scale as < > ~ 1 2 for long times, 
i.e. sub diffusively. This behaviour corresponds to the naive dynamical expo- 
nent z = 4 in model B (in the nomenclature of Halperin and Hohenberg p3| ) 
and has been expected, for the system here is critical with respect to the 
transverse directions and no renormalizations are necessary in this subspace. 
The extrapolation of the anomaly-exponent 77 into low dimensions is difficult 
in view of the high upper critical dimension dc = 9. Naively, we can simply 
set e = dc — d = 6 in Equation (^) to predict rj i.e. in a 3-dimensional 
system, resulting in an estimate 

ri(d = 3) = 1.20 . (44) 

Although e has been considered as small quantity the two loop-correction 
with respect to the one-loop result is merely 20%, even at e = 6. Due 
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to these small corrections we expect that equation (^), even though it is 
strictly valid only near dc = 9, produces good approximations for t] also in 
low dimensions. 



3.2 Random Potential With Equally High Potential 
Mountains 

This model is analyzed in analogy to the model with unsymmetric random 
potential of the last section. The quenched disorder according to Equation 
([l3| ) resulting from the random potential with equally high potential moun- 
tains leads, in analogy to Equation (]T3|), to the dynamic functional 

J[s,s]=Jd'^r S^Jdt[s (s + A(A^(kxA± - +pA||(K||A|j - 



1 2 



A dtA±s 



— a 



A dtVws 




that still contains irrelevant terms to be eliminated. 

Transverse and longitudinal lengths scale as before r± ~ ry ~ . 

Since the dynamic functional is dimensionless the dimension of fields and 
coupling constants are 

-4 



At 



d-S 

7-d 
fl 2 



d+5 
/i 2 





a ^ fi" (7 ~ yU 2 . (46) 

Thus, the coupling constants k,\\ and n are again irrelevant in the renormal- 
ization group sense, but a is not. This signifies that in the given random 
potential both longitudinal and transverse noise are relevant. The dimension 
of the nonlinear coupling g shows that 

dc = 7 (47) 

is the upper critical dimension of this model. The finite ry and the transverse 
coupling constants a and k± are again absorbed by a suitable scale change. 
Thus, the appropriate dynamic functional for the transverse phase transition 
in the given random potential is 



d'^r 



dt 



ss + Xs{A±{A± - T±] 


- pA 


l)s + ^A(7(V| 






2 








X JdtA^s 


— a 


A JdtViis 


■] 



PT^ „2 



(48) 
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This dynamic functional exhibits the same symmetries as the one in the case 
of the unsymmetric random potential, except the longitudinal scale transfor- 
mation, where here the parameter o is additionally transformed according 
to cr — > In addition to the parameter combination g^p~^ ^ ^P~^ is also 
an invariant variable of this model and the perturbation expansion will be 
organized in powers of both variables. 

Due to the quenched disorder, the Gaussian correlator is again time inde- 
pendent so that, for simplicity, we also here transform J into a quasi-static 
Hamiltonian via fE6D and (|2T 



+(^(-A^ + aA|| 



In comparison with the former model, the vertex and the Gaussian propaga- 
tor remain the same (pQ]), whereas the Gaussian correlator of the quasistatic 
model reads 

" |qi(qi + ri)+p«|p- ^ ' 

By dimensional analysis taking into account the causality properties we find 
the vertex functions Fi^i, r2,o, and Fi 2 primitively divergent. 
Neglecting higher orders in q and e = dc~ d, a one-loop calculation for these 
vertex functions, with dimensional regularization, results in 

fi,i(q) = qi(qi + rx) + pgj + ^$gjr:^ (51) 

e p2 

1 B 

f2,o(q) = -2ql - 2(TgJ - t — — gf^i 

4 e p2 

ri,2(q) = ^9(l\\-^-: ——(l\\^±' 
4 e p2 

Be = T^Od-iT(l + -)T(^^)T(-) (52) 

is a suitably chosen e-dependent factor. Notice that although not performed 
a two-loop calculation could here also be done with the method described in 





where 
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the Appendix. 

The primitive divergences are absorbed in the redefinition of the parameters 



p = Zpp, a = Z„a, g = p^Z^u 



(53) 



In comparison to the latter model the coupling constant a has additionally 
to be renormalized. From Equation (^) we easily obtain the Z-Factors in 
minimal subtraction 



1 



V 



t 

\v 1 



(54) 



8 e w 

l + l^(l + ^) + 0(t;2), 



expressed as functions of the dimensionless renormalized parameters v := 
B^u^p~2 and w := ap^^ that are invariant under the longitudinal scale trans- 
formation. 

In analogy to the former model we obtain the renormalization group equation 



^ d „ (9 . d d 

+ Pw^ + PC^ + P^ 

ov ow op op 



rn,n({g|i,q±},r±,t;,w,p,/i) = (55) 



with the parameter functions depending on v and w 



Pv 



dv 
dp 



dw 

P. := P^ 



p- 



dlnp 



dp 



'-v{4: + w) + 0{v^ 



--v[5 -6w + w^] + 0{v' 



(56) 



The associated characteristics are defined by 



dl 



lj^lnpil)=Cm,wil)) (57) 



with the initial condition w{l = 1) = w and the initial conditions from 
Equation (|36|) . In the scaling limit / <C 1, v{l) and w{l) fiow to an infrared 
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stable fixed point {v^,w^) given by the zeroes of and with a positive 
gradient. The zeroes of (3^) in this order are 

Wi = 1 W2 = 5 , (58) 

and those of /3„ 

v, = -2 = jf^ • (59) 

This yields the infrared stable fixed point 

w, = l + 0(e) v, = '^e + O(e^) . (60) 



The domain of attraction of this fixed point can easily be recognized as 

v>0 w <5 + 0(e) . (61) 

In the case w > 5 + 0(e) the critical behaviour of the system is dominated 
by the degenerate fixed point 

= oo = . (62) 

Concerning this degenerate fixed point we remark the following: 

i) For w ^ oo the transverse noise vanishes in comparison to the longitu- 
dinal one. This can be seen easily by substituting pgy gy and extracting 
ap~^ = w from the correlator. Then the remaining coefficient of the lon- 
gitudinal part in the numerator of the correlator is 1, the coefficient of the 
transverse part is w~^. At the degenerate fixed point, this system here be- 
haves as the model with a random potential with equally deep potential 
valleys and whose noise is therefore given by Equation (0). This model is 
analyzed in the next section, but we anticipate some results concerning the 
degenerate fixed point. According to Eq. (|76D the degenerate fixed point 
possesses a finite fixed point value {v ■ w)* = |e + C(e^). Although the de- 
viations from normal diffusive behaviour only appear in order two-loop, the 
positive r] from Eq. ( [77D demonstrates the system to be also superdiffusive 
at the degenerate fixed point. 

ii) Both a "normal" and a degenerate fixed point are also observed in the 
noncritical model with quenched disorder [0]. While here at the transverse 
phase transition the situation of two fixed points appears in the random 
potential with equally high mountains and the degenerate fixed point is de- 
scribed by the model with the random potential with equally deep valleys. 



19 



in the noncritical region, however, the situation of a normal and a degener- 
ate fixed point occurs in the model with unsymmetric random potential and 
the degenerate fixed point is described by both the model with equally high 
mountains and equally deep valleys, for these latter models are identical in 
the noncritical region. 

Now we proceed to investigate the normal fixed point (|60D. For this fixed 
point the anomaly-exponent reads 

V = l^ + Oie') . (63) 

At the fixed point, the second characteristic has the form 

p{l) = p/-2. (64) 

in analogy to Equation (^Of) . 

Returning to the dynamic model and again exploiting the renormalization 
group equation at the fixed point, dimensional analysis, and the invariant 
scale transformation we obtain the universal scaling behaviour of the vertex 
functions in the asymptotic limit 

rn.nllgihqx,^^},^^,^^*,^^*, A,p,/i) = (65) 

^-i„{n+n-2)-iri(d+5)-ln(d-3)+d+5p_ / JL Si ^\ 1^ y yj X n n] 

Thus, only longitudinal lengths scale anomalously below the upper critical 
dimension dc = 7. From this equation we derive the scaling form of the 
density response function 

X(q,t) = f{q!t'^'^',qlt'^) (66) 

which coincides with the form found for the preceding model. The positiv- 
ity of 1] = |e + C'(e^) signifies superdiffusive behaviour below dc = 7. A 
comparison with the corresponding model without quenched disorder again 
shows the quenched disorder to be the reason for the enhanced spread of 
fiuctuations in the driving force direction in a dimensional interval that is 
here, however, 5 < d < 7. 

A straightforward extrapolation of the anomaly-exponent into three dimen- 
sions by setting e = 4 into the one-loop result gives the approximate value 

r]{d = 3) = ^ . (67) 
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This numerical value is very distinct from the one-loop and two-loop values 
for T) in the model with unsymmetric random potential. 
We now discuss the case that the transverse noise term with 7 7^ of Equa- 
tion ([T5|) is generated by coarse graining even in this symmetric random po- 
tential. Then an additional scaling variable j/l'^ arises in (^). The crossover 
exponent is here = 2 and it simply reflects the naive dimension of the cou- 
pling constant 7, because 7 itself is invariant under the longitudinal scale 
transformation and for dimensional reasons the transverse noise term with 
7 7^ needs no additional renormalization. Since 7 is a relevant variable, 
for 7 7^ the system eventually flows to the fixed point of the model with 
unsymmetric random potential. 



3.3 Random Potential With Equally Deep Potential 
Valleys 

Having no transverse noise according to Equation ( p!4D this model is formally 
obtained from the model with random potential with equally high mountains 
by setting the transverse noise coefficient a = 0. Then the relevant dynamic 
functional is here 



d\\ dt 



1 



ss + As(A_l(A_l - Ti_) - pA||)s + -A5((V||s)s^ 



—a 



A dtVns 



The dimensions of lengths, fields, and coupling constants as well as the upper 
critical dimension dc = 7 are as in the latter model. The symmetry properties 
of this dynamic functional are the same as in both preceding models with the 
exception of the invariant scale transformation. The lack of the transverse 
noise term has the consequence that the dynamic functional is invariant under 
a scale transformation depending on two parameters 



/3r|| 



as 
2 , 



r± 
s - 

9 - 



a-^(3-^s 
ap'^g 



(69) 



a ^/3(T . 



Therefore, the combination of coupling constants g^ap 2 is the appropriate 
invariant variable in this model and is exactly the product of the variables 
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being invariant each in the latter model, but not here. 

Because of the quenched disorder we again use the transformation to a quasi- 
static Hamiltonian 

n[ip,^]= y'rf^r|(^(Ax(Ax-rx)-pA||)(^ + ^(7(V|i<^)<^' + a(^A||(^| . 

(70) 

Vertex and Gaussian propagator are the same as in both preceding models, 
whereas the Gaussian correlator reads 

2aql 

^ [qi(qi + Tx) + pgj}]2 ' ^^^^ 

As in the previous model Fi i, r2,o, and Fi 2 are primitively divergent. We 
have computed them in the lowest nonvanishing order, i.e. we have per- 
formed a one-loop calculation for r2,o and Fi 2 and a two-loop calculation 
for Fi 1. The two-loop calculation involves similar integrals as in the model 
with unsymmetric random potential and has also been performed with the 
help of inverse Mellin transformation (cf. Appendix). 

In dimensional regularization we obtain the primitively divergent vertex func- 
tions, up to higher order terms in q and e = dc — d 

fi,i(q) = qi(qi + r^)+/3gJ-4:^^g2^-(0 - 0.001799e) (72) 



P 

2 J- -^e y 2 



r2,o(q) = -2ag|| - - — ^^g||r^ 



4 e 



ri,2(q) = ^9q\\ — —q\\r±^ , 
4 e p2 

where -B^ is defined by Equation (|52D . The same redefinition of the coupling 
constants as in Equation (^) absorbs these divergences and yields in minimal 
subtraction 

Zp = 1 - 0.007196-^ + 0(1;^) (73) 
o e 

Zu = l + ]- + 0{v'), 
4 e 
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where v := B^u^ap~^ is the dimensionless renormahzed variable of the modeL 
This variable is invariant under the longitudinal scale transformation (Eq. 
(p9|)) and is the product of v and w in the latter model, but it will here 
merely be denoted by v for simplicity. 
The renormalization group equation reads 



d_ 

' dv 



d 



d_ 



d_ 
dp 



rn,n({g||,q±},rj_,t;,p,a,/i) = (74) 



with the Wilson parameter functions 

dv 



Pv 



o 

din a 



dp 
dlnp 



dp 



o 

-0.0U39v^ + 0{v^ 



(75) 



only depending on v. 

The infrared stable fixed point, being a zero of is found at 



8 



Hence we obtain the anomaly-exponent 

7] = 0.0512e2 + 0{e^) 



(76) 



(77) 



The characteristics of the renormalization group equation are defined by 
Equation (^5|) and an additional equation for a{l) that has the same structure 
as the equation for p{l). At the fixed point the solutions for these flow 
equations read 

p(/) = p/-2^ a{l) = af^* , (78) 



where 



(79) 



Returning to the dynamic model and combining the renormalization group 
equation at the fixed point, dimensional analysis, and the invariant longitu- 
dinal scale transformation we obtain the universal scaling behaviour of the 
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vertex functions in the asymptotic limit 

rri,n({gii,q±,Lu},r^,^^*, A,p,cr,/i) = 

^-r;(-l+|n-ifi)+lC<T*("-")-|"(rf+5)-5n(d-3)+d+5 

We again see that only longitudinal lengths scale anomalously. In contrast 
to both preceding models there is here an additional exponent appearing 
in the global /-factor. As this exponent is cancelled for Fi^i, the density 
response function nevertheless has the scaling form 

X(q,t) = (81) 

in analogy to both preceding models. As the anomaly-exponent r] = 0.0512e^+ 
(!?(e^) is positive and the upper critical dimension rfc = 7 is the same as in the 
preceding model, all statements concerning superdifFusion and speeding-up 
of fluctuations by quenched disorder are here also valid. 
The straightforward extrapolation of rj from dc = 7 to d = 3 hj setting e = 4 
produces the approximative value 

ri{d = 3) = 0.82 . (82) 

This numerical value is close to the one-loop value for 77 in d = 3 in the 
model with random potential with equally high mountains, but is very far 
from the one-loop and two-loop values in the model with unsymmetric ran- 
dom potential. 

We now investigate the possibility that even in this symmetric random po- 
tential transverse noise terms are produced by coarse graining. First, we 
consider the case that transverse noise proportional to a / d'^r[\ J dtA±s]'^ is 
generated. Then we are in the situation of the preceding model. If a is small 
the system is in the region of attraction of the degenerate fixed point dis- 
cussed in this section. This implies that the transverse noise vanishes under 
the renormalization fiow and the results of this section remain valid. 
Second, if transverse noise proportional to 7 / d'^r[X J dt'V±s]'^ is produced, 
an additional scaling variable j/l^ arises in (^). Due to dimensional rea- 
sons, this relevant operator needs no additional renormalization. The invari- 
ant scale transformation according to Equation implies 7 — > /5~^q;~^7. 
Thus, the crossover exponent is related to the exponents r] and Co-* by = 
2 + 2?7 — For 7 7^ the system eventually fiows to the fixed point of the 
model with unsymmetric random potential, because 7 is a relevant variable. 
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4 Longitudinal Phase Transition 



The three different reahzations of the random potential I-III described in 
Section 2 converge to a single model in the region of the longitudinal phase 
transition as is sketched in the following. 

For the three realizations of the random potential the dynamic functional 
still containing many irrelevant terms is given by Eq. ([T5|) and by Eq. (^51) 
with a 7^ and a = 0, respectively. The longitudinal phase transition is 
characterized by ry — > and finite t±. The external scale fi"^ here measures 
small ry. Comparing the leading gradient terms in driving force direction 
(~ AspA||K|| A||s) and transverse direction (~ Asr_|_A_|_s) we find for length 
scales 

r\\ ~ fi^^ r± ~ . (83) 

As the dynamic functional is dimensionless the coupling constants of the 
transverse noise scale as 7 ~ a ~ fi~^, i.e., only the longitudinal noise is 
relevant. Hence, only a single model is required to describe the longitudinal 
phase transition, independent of the kind of the random potential. Further, 
the dimensions of fields and the other coupling constants are 



At ~ /i~ 



5( ~ /i 2 - 



(84) 



indicating that k,± and k are irrelevant and that 

dr = 6.5 



(85) 



is the upper critical dimension of this model. 

After a suitable scale change of fields and lengths the relevant dynamic func- 
tional for the longitudinal phase transition is 



J[s, 5] 



rfV dt 



ss + As(-A_L + pA||(pA 



\\-rii))s + -Xg{V\\s)s^ 



-a 



A dtVws 



(86) 



The symmetry properties are the same as in the last model of Section 3 
inclusively the invariant longitudinal scale transformation, because transverse 
noise exists in neither model. As a consequence g'^(7p~2 is the appropriate 
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invariant variable of this model. 

Due to the quenched disorder, we transform to a quasistatic Hamiltonian 

n[ip, = jd'^r {<^(-A^ + pA|| (pA|| - r||))99 + \g{V^\^)^^ + a<^A||^] 

(87) 

Prom this equation we directly read off the elements of perturbation theory, 
i.e., the vertex —igq\\ and the Gaussian propagator and correlator 

^ qi + pgj(pg| + r||) ^^^^ 

[qi + pgf(pgj + r||)]2 " 

The vertex functions Fi i, r2,o- and ri^2 are primitively divergent. The sin- 
gular parts of Fi 1 are proportional to and the singular parts of r2,o 
and Fi 2 are proportional to and gy, respectively. 

In a one-loop calculation we obtain the primitively divergent regularized 
vertex functions in bare quantities up to higher orders in q and e — dc — d 

fi,i(q) = ql + pgj(pgj + f||) + -^|^gjrf (-pgJ + 0-f||) (89) 

f2,o(q) = -2agJ - i^^,Jff 

D e p2 " 

ri,2(q) = m\\-^-^——(i\\rii , 

where 

= (^^^-ir(i + e)r(il^)r(^) (90) 

is a suitably chosen e-dependent factor. These divergences are absorbed in 
the redefinition of the coupling constants 

p = Zpp f|| = Z^||T|| & = Z^a g = ii'^ZuU . (91) 

In minimal subtraction we obtain the Z-factors as a function of the dimen- 
sionless renormalized invariant variable v :— C^u^ap~2 

Zp = + (92) 
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Zu 



1 + i- + 0(v'-) 

D e 



Here, the renormalization group equation is of the form 

d ^ d ^ d d 



ov op oa ary 



d_ 
dp 



^^^,„({g||,q±},r||,^;,p,c^,/i) = 

(93) 



with the Wilson functions 



dv 
dp 
dhia 



-2v 



e v + Oiv^ 

24 ^ 



(94) 



dp 
dlxip 



dp 
d In T\\ 





b 



dp 



6 



t; + C(w2) 



In addition to the flow equations for v{l), p{l), o"(/), and p{l) being of the 
same form as in the preceding model there is a flow equation for f||(/) that 
reads 

|lnf,|(/) = «:(^;(/)). (95) 
In the scaling limit / ^ 1 that corresponds to the relations 



< 1 



< 1 



p^ 



< 1 



(96) 



v{l) flows to an infrared stable fixed point which is directly obtained as a 
zero of Py 

'^^ 2 
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V* = —e + U{e-} . (97) 

At the fixed point the solutions of the flow equations are given by Equation 
(78) and by f||(/) = r||/'** with /t* := k,{v^:). Inserting the fixed point value 
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for f * into (p we find the anomaly-exponent 

V=^e + 0{e') . (98) 

The fixed point values of the other parameter functions are 

U = -^e + Oie') «:, = le + 0(e2). (99) 

Returning to the dynamic model and combining the solutions of the renor- 
malization group equation at the fixed point with a dimensional analysis and 
the invariant scale transformation we finally arrive at the universal scaling 
behaviour of the vertex functions at the longitudinal phase transition 

rn,n({g||, q±, r^, v^, A, p, a, ij) = 

^-»?(-l+|n-in)+|C<T.(n-n)-5n(2a!+5)-|n(2d-7)+2d+3 

/ f gil q I CiJ 1 Til \ 
■Tn,n [[jTT^^ "F' /4 / ' ^4-' ^' P' ^' /^j (100) 

In addition to longitudinal lengths the critical parameter ry also exhibits 
anomalous scaling behaviour. In analogy to the models in Section 3 we 
derive the scaling form of the density response function 

X(q,t) = f{qp^^^'^'\qlt) ■ (101) 

While the system is normally diffusive with respect to the transverse direc- 
tions, in the critical longitudinal direction typical length squares scale for 
long times as 

< rj > ~ , (102) 

i.e. in comparison to the naively critical (corresponding to the naive dy- 
namical critical exponent z = 4 from model B ||23|) the spread of fluctuations 
in the driving force direction is enhanced, due to the positive r] = ^e + (9(e^). 
A comparison with the corresponding model for the longitudinal phase tran- 
sition in a system without quenched disorder p shows the surprising result 
that the longitudinal phase transition is here continuous, evidenced by the 
existence of an infrared stable flxed point, whereas in the model without 
quenched disorder there is no infrared stable flxed point. 
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5 Summary and Outlook 



We have analyzed the transverse and longitudinal phase transition in uni- 
formly driven diffusive systems with quenched disorder. These systems show 
a wide variety of possible scenarios, because the symmetry properties of the 
random potential are an additional distinguishing feature to which universal- 
ity class a model belongs. In the region of the transverse phase transition the 
three different random potentials I-III (Fig. |^) actually define three different 
models and in the noncritical region they still define two different models all 
of which are part of different universality classes. 

Together with earlier investigations the renormalization group studies of an 
entire model class are hereby completed. This model class includes the driven 
diffusive systems with and without quenched disorder both in the critical re- 
gions of the transverse and longitudinal phase transition and in the noncrit- 
ical region. Fig. ^ gives a graphical overview of the model class, where the 
single models are ordered chronologically from left to the right. The models 
in an ordered substrate i.e. without quenched disorder, had been stud- 

ied prior to this work for all three regions of the phase diagram mentioned 
above. The models for a system with quenched disorder in the noncritical 
region |^ had also been investigated. The other models of Fig. ^ have been 
studied in the present paper. 

The upper critical dimension dc is different from model to model and varies 
from 2 to 9. Below dc-, the vertex functions being typical statistical quantities 
of such systems show universal anomalous scaling behaviour on large length 
and time scales. The deviation from pure diffusive behaviour is characterized 
by the anomaly-exponent rj. It indicates how strongly longitudinal lengths 
scale anomalously. In Fig. |^ the result for r] in the highest calculated loop 
order is given for every model. Note that, due to Galilean invariance, rj is 
even exact in all models without quenched disorder. 
We emphasize the following results: 

In all models of this model class (except for the longitudinal phase transition 
in a system without quenched disorder) the anomaly-exponent rj is positive 
implying superdiffusive spreading of density fluctuations in the driving force 
direction. 

A model with quenched disorder always has a higher upper critical dimension 
than the corresponding model without quenched disorder. Due to the field 
theoretic results, in the dimension interval between these two upper critical 
dimensions the quenched disorder is the reason for superdiffusive spread of 
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density fluctuations in the high temperature region and at the transverse 
phase transition, respectively. Notice that the anomalous diffusion of fluc- 
tuations is not directly connected with the behaviour of transport coeffi- 
cients relating the mean current to the mean density, because a fluctuation- 
dissipation theorem (Einstein relation) does not hold in this strong nonequi- 
librium situation with quenched disorder. Within this model class we have 
not found a qualitative argument why disorder generates super diffusion. We 
mention, however, a study of a one dimensional driven lattice gas p^,|25[] 
where quenched disorder is not, as here, spatially fixed, but associated with 
moving particles. The authors found a superdiffusive spread of a jam behind 
the slowest particle and a superdiffusive spread of free spacing in front of 
it. Whether this situation can be transferred to our models by identifying 
the deepest potential valley and the highest potential mountain, respectively, 
with the slowest particle, should be a topic of further investigations, espe- 
cially Monte Carlo simulations. 

For the longitudinal phase transition there is no uniform statement. In the 
model without quenched disorder no infrared stable fixed point has been 
found and some analytic arguments point to a discontinuous phase transition, 
whereas an according two dimensional lattice gas demonstrates a continuous 



phase transition in Monte Carlo simulations p6[. In the model with quenched 



disorder studied here, however, we find an infrared stable fixed point and thus 
a continuous longitudinal phase transition. 

Despite the partially high upper critical dimensions it appears that the 
anomaly-exponent (when existing) may be extrapolated quite accurately into 
low dimensions. Two observations lend support to this procedure: first, the 
two-loop correction to 1] in the critical transverse model from Section 3.1 is 
small and second, the coefficients of 1] are small in all models. 
For the two-loop calculation in critical models with quenched disorder a 
new technique has been developed enabling us to manipulate mixed q^- 
propagators and -correlators. This technique is based on an inverse Mellin 
transformation and is described for one model at the transverse phase tran- 
sition in the Appendix. This method is applicable to all critical models with 
quenched disorder of this model class, but has only been performed for two 
models (Sec. 3.1 and 3.3). Furthermore we expect it to be useful for two-loop 
calculations in other physical problems where q^-propagators are involved. 
Monte Carlo simulations of two-dimensional driven diffusive lattice gases 
without quenched disorder ||I],^[Tl|,|12[ are in excellent agreement with field 
theoretic predictions for the transverse phase transition and the noncritical 
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region [^^. For the transverse phase transition in systems with quenched 
disorder, however, there is a simulation study ||2^ that is, for two reasons, 
hardly compatible with the models investigated here by field theory. First, 
the quenched disorder was there modelled by randomly blocked sites and not 
by random potential barriers between the sites. Second, the concentration of 
blocked sites is so small that only the crossover behaviour between a system 
with and without quenched disorder was observed. Further, we mention a 



recent simulation for a noncritical system with quenched disorder |28 



It is desirable to compare the field theoretic results obtained here for the 
various critical models with corresponding Monte Carlo simulations that are 
still to be done. These simulations are a nontrivial challenge, because besides 
the average over a huge number of realizations of the quenched disorder it 
is to pay attention to the fact that the periodic boundary conditions usu- 
ally used let pass a particle repeatedly through the system and let it see the 
same quenched disorder as before. By this unwanted correlations of the ran- 
domness enter into the simulational results that make it difficult to compare 
them with the field theoretic results that are based on the assumption of 
uncorrelated disorder. Therefore Monte Carlo simulations of driven diffusive 
systems with open boundaries already done for systems without quenched 
disorder p9[ seem to be more appropriate. 

We finally remark that we have extended the model class investigated here 
in the way that we allow for random particle sources and drains in the diffu- 
sive systems. The noncritical model without quenched disorder with such a 
particle nonconserving randomness had already been studied More- 



over we have analyzed the influence of nonconserving noise onto the critical 
behaviour in systems with and without quenched disorder which is demon- 
strated in a paper soon to be published fl^ . 
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A Two— Loop Calculation of Fi i 



In the Appendix we explicitly show the two-loop calculation for the model 
that describes the transverse phase transition in an unsymmetric random 
potential (cf. Section 3.1). The graphical elements of the diagrammatic per- 
turbation expansion are shown in Fig. § directly for the quasi-static model 
which is here used to facilitate the calculation. The cp-legs are indicated by 
an arrow and the q\\ of the vertex by a dash perpendicular to the propagator 
line. The mathematical expressions for the graphical elements are given by 
Equation (H). 

In this model we have to calculate the primitively divergent vertex functions 
Fi 1 and Fi 2. There are eight two-loop diagrams contributing to Ti i (Fig. 
^ and 25 two-loop diagrams contributing to Fi^2 each of which obeys causal- 
ity that forbids closed propagator loops. Momentum conservation demands 
that at each vertex the sum over all wave vectors is zero. Evaluation of these 
diagrams requires integration over all internal wave vectors. 
While the one-loop calculation is easy to perform by standard methods, 
the two-loop calculation involving mixed q'^-propagators requires more so- 
phisticated tools and has not been described in literature. We introduce a 
technique that is based on an inverse Mellin transformation. 
For Fi 1 we have to compute the diagrams b[^'^^ to Bg''^^ from Fig. |^. We de- 
note the external momentum by q and the internal ones to be integrated over 
by p and k. As Fi^i is quadratically divergent because of CP-symmetry and 
dimensional reasons and the external (p already provides a factor q\\, the parts 
of the integrands being proportional to q\\ contain all singularities. Therefore, 
the integrands are first expanded with respect to the external momentum gy 
to first order. For simplicity, we now substitute p^p\\ — * p\\, p^^y — >■ k\\ and 
from now on we omit the superscript " ° " to characterize unrenormalized 
quantities. 

In the next step all even powers of p\\ and fcy in the numerator of the inte- 
grands are written as p| = [pi(pi + T±) — pi(pi + T±) {kj analogously). 
The first term of the right hand side cancels against factors in the denomi- 
nator. For all odd powers the half of the integrand is reflected with respect 
to k\\. This step reduces the superficial degree of divergence of the p- and 
k-integration by 2. It is allowed because the integration runs over the whole 
/cii-axis. 
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Then all integrals are reduced to the two types of integrals 

I{a,(3,-f,6,^i,iy) := (103) 
Pi" f (Px-kx)^^ ki^ 



P [Pi(pi + rx) +Pir [{P± - k±)^ + {Pii - [ki + k 



2]u 



f P i Pi 

Jp [Pi(pi + 

hkr r (p±-k^)27 {p±-k^r 



k [ki + kj 



[(p± - + (Pli - W [iP± - k±)^ + {P\i + h 



l21<5 



the arguments of which can be 0, 1, 2,... 

As the integrands are already expanded with respect to gy, the whole inte- 
gration over p and k is at most logarithmically divergent. The integration 
variables are chosen such that the p-subintegration is always primitively con- 
vergent, i.e., the naive dimension 5-p of this integration (measured in powers 
of the external momentum scale /i) is negative. In the integrals of both types 
the k-integration is at most logarithmically divergent. 

Due to these naive dimensions of the p-, k- and the whole integration it is 
possible (and necessary for the following calculus) to set tj_ = in the k- 

integration, because the r_|_ 7^ 0-parts only provide convergent contributions. 
The sum of the two-loop diagrams for Fi 1 expressed by the integral types / 
and F reads 

8 4 

^Sp) = 4^gJ[3/(l, 3, 0,1, 1,2) -3/(1, 3, 0,1, 3, 3) -2/(3, 4, 0,1, 1,2) 
i=i P 

+ 2/(3, 4, 0, 1, 3, 3) - 2t^/(2, 4, 0, 1, 1, 2) + 2t^/(2, 4, 0, 1, 3, 3) 
+ ^/(O, 2, 1, 2, 1, 2) - ^/(2, 3, 1, 2, 1, 2) + /(4, 4, 1, 2, 1, 2) 

- F(3; 1, 4, 0, 1, 1, 3) - \f{\- 1, 3, 0, 1, 1, 3) + F(l; 1, 4, 1, 2, 0, 1)]. 

(104) 

By factorizing the q*-denominators in / und F into transverse and longitu- 
dinal parts we are in the position to apply the successful methods that are 

used for q^-propagators and -correlators. This factorization is done by the 
Mellin transformation. We demonstrate the method for the integral type / 



33 



in detail. The calculation of F goes analogously. 

The Mellin transformation of the function (a + is 



oo 



dxx'-\a + x)-" = r(^)r(a t) ^,_^ ^ ^^^^^ 



where the conditions a > and < Re{t) < Re{a) must be fulfilled ||3T 
The corresponding inverse Mellin transformation 

where the integration path parallel to the imaginary axis is restricted by 
< to < -Re(a) [0, proves to be the appropriate tool to factorize the 
denominators in /. 

First we only treat the k-integration of the integral type / (|103|) and apply 
the inverse Mellin transformation to both denominators 



k[(p±-k^)4 + (p||-fc|l)2]5[ki + A;f]- 
to+ioo dt p+ioo T{t)T{6 - t)T{s)T{iy - s 



to-ioo 27rz Jso-ioo 2m T{5)T{iy) 
1 1 



i 



k [fcf]1(j9|| - A;||)2]* [ki]2(-^)-A'[(p^ - k^)2]2{5-t)-7 " 

In the integrand transverse and longitudinal momenta are now separated and 
only quadratic. The k-integration can now be performed with the help of 
the usual Feynman relations 

1 1 nOO 

dss'^-^e-'^ (108) 



A" Via) Jo 



1 _ r(E.«.) r^idx,^?^-' ^^^^^^-1) 



where T{a) is Euler's F function. The result reads 

^ 2{2ttY ^2' ^ 2 ' Jto-ioo 2mJso-ioo 2m r(5)F(z/) 

T{s + t- |)F(| - g)F(| - t) T{2u + 2(5 - /i - 7 - ^ - 2(g + t)) 
F(s)F(t)F(l - s -t) F(2(z/- s) -/i)F(2(5-t) -7) 
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r(^ - 2{v -s)+ /x)r(^ - 2{5 - t) + 7) 
r(d - 1 - 2z/ - 25 + /i + 7 + 2(s + t)) 

The k-integral exists under the conditions 

2(so + to) > 1 (110) 
1 1 
«o<2 < - 

4(z/ + 5- (so + to)) -2(yU + 7) > 

- 1 - 4(z/ - So) + 2/i > 

rf- 1 -4((5-to) + 27 > 

which prevent UV and IR divergences, respectively, at the k\\- and k_L- 
integration. These conditions restrict the complex integration paths of the 
s- and t-integration which are here dependent on the spatial dimension, due 
to the dimensional regularization. 

The result of the k-integration is according to Equation ( |109| ) proportional 
to powers of and |p_l| with exponents that are dependent on s, t, and d. 
Thus, the remaining p-integration is analogous to the one-loop problem, up 
to changed exponents. It is straightforwardly performed and we obtain for / 

r(2(/j + 5 + z/)-(a + 7 + /x)-((i+l)) 

I{a,P,-f,d,fX,u} = Cd 2(/3+<5+.)~(a+7+A^)-(rf+l) 

•/ 7^ — r(2z/ + 25-^-7-^-2(s + t))r(s + t--) 



2 ' '2 ' r(rf-l-2z/-25 + /i + 7 + 2(s + t)) 

fllll 



r(^ - 2[v - s) + /x)r(^ - 2(5 - t) + 7) 



r(2(z/-s)-/x)r(2(5-t)-7) 

-, fQ \2 r(i)r(^— i-) 

where Cd = ^ (12^) v{p)v{8)v{v) ^ constant that is different for every / 
and depends on the dimension. The conditions for the existence of the p- 
integration are 

l-(so + to) > (112) 
c? + so + to + a + 7 + /i-/5-2(5-2z/ > 
2(/3 + 5 + z/)-(a + 7 + /i)-(rf+l) > 0. 
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After the momentum integrations there remain parameter integrations with 
respect to the Melhn variables s and t over parallels to the imaginary axis. 
The complex integration paths and the spatial dimension d are restricted by 
the conditions (|110| ) and (|112|) . With e = 9 — d we obtain from the first four 
inequalities of (|110|) 

^<So + to<'5 + z/-^(7 + ^)-2 + |, (113) 

where sq < |, to < |- All other conditions are satisfied for all / of Equation 
([10^) , if < e < 2. 

Since the sum sq + to appears in the inequality (|113|) we transform to the 
new integration variables 

z:=s + t w:=^{s-t). (114) 

The integration paths are now determined by 

^ < zo < 5 + z/-^(7 + /i)-2 + | (115) 

kol < \. 

After this substitution Equation ( [111| ) gives 

T( n , ^ r. r(2(/3 + ^ + z/)-(a + 7 + /x)-(rf+l)) 

rzo+ico rwo+ioo dyj ^ d—1 „ 1, 

' ZQ—ioo 2,7ri Jwo—ioo 

where the abbreviation 



/ 77- / Tr-r(2^^ + 25 - /X - 7 - —- - 2z)T{z - -)f{z, w; e) 



f{z,w;e) := r(5 - ^ + ^)r(z. - 1 - ^)r(- - 1 - t.)r(- - ^ + ^) 

r(^ -2{u-^-w)+ fi)T{^ - 2((5 - f + w) + 7) 

r(2(z/-|-«;)-^)r(2((5-f + «;)-7) 
r(9 - e + z + a + 7 + /i - /3 - 25 - 2z/) 



r(8 - e - 2z/ - 25 + + 7 + 2z) 

denotes the part of the integrand that is free of poles. 

To extract the divergent parts of / we have to distinguish two cases. 



(117) 
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Case 1: The k-integration is primitively convergent. 
In this case the parameters satisfy 

i/ + (5-^(7 + /i) = 3, (118) 

which is true for /(O, 2, 1, 2, 1, 2), J(2, 3, 1, 2, 1, 2), and /(4, 4, 1, 2, 1, 2) from 
Equation ( p.04] ). According to Equation ( p.l5| ) this means for the constant zq 
which determines the integration path 

l<zo<l + l. (119) 

For these / the whole integration over p und k is logarithmically divergent 
so that their parameters fulfill the equation 

2(/5 + 5 + z/)-(a + 7 + /i)-(d+l) = e. (120) 

Thus, the coefficient r(e) = -r{l + e) of the double integral over z und w 



(|116|) contains an e-pole, whereas the double integral itself is convergent, as 
the integrand is free of poles even in the limit e — > and the real part of the 
integration path can be chosen between | and 1 according to Equation (|119|) . 



The double integral in Equation ( p.ltj| ) can therefore be computed at e = 0, 
because together with the e-pole as coefficient only convergent contributions 
are produced for e 7^ 0. The convergence of the double integral is ensured by 
the asymptotic behaviour of the F function with complex arguments [p2| 



lim |F(a; + %)| e^l^l = {2n)^ . (121) 

Ij/I^oo 

Due to the accumulation of F functions we are not able to perform the inte- 
grations with respect to z and w analytically. Hence, each of the three / is 
computed numerically for e = 0. 

Case 2: The k-integration is logarithmically divergent. 
In this case the parameters are restricted to 

z/ + 5-i(7 + /i) = ^, (122) 

which is correct for the remaining I from Equation (|104|) . According to 
Equation ( |115| ) the integration path is here restricted by 



i<--o<i + j. (123) 
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In the limit e ^ the integration path is trapped by this condition between 
the poles of the integrand ( |116|) at z = ^ and -2^ = | + f • In order to extract 
the e-poles the complex integration path is decomposed into two parts (as 
shown graphically in Fig. ^ : The first part is a line parallel to the imaginary 
axis with an arc to the left of the singularity z = ^, while the second part is 
a circle around z = ^. The z-integral over the circle gives the residuum at 

This way we obtain from Equation ( p.l6| ) 

a , ^ ^T{2{P + S + u)-{a + -f + fi)-{d+l)) 
I{a,P,j,6,^i,u) = C,^^ 2(/^+i+.)-(a+,+,)-(Zi) (124) 

V{1 + - - 2z)V{z - -)f{z,w-e) 



ioo 2712 Jwo~ioo 2,71% 2 2 

e 



WO+ioO don \ 

+ r(-) / —f{z = -,w-e) 

Naturally, the integration path of the ^-integration is again chosen as a par- 
allel to the imaginary axis, whose real part Zq now lies between — | and |. 
After the integration path is changed, both the double integral over z and 
w and the single integral over w are eventually convergent. Only the coef- 
ficients of the integrals contain the e-poles. For the further calculation we 
have again to distinguish two cases, 
a) The whole integration is logarithmically divergent. 

This is the case for the integrals /(I, 3, 0, 1, 1, 2), J(l, 3, 0, 1, 3, 3), /(3, 4, 0, 1, 1, 2), 
and J(3, 4, 0, 1, 3, 3) from Equation( |l(J4 ). Their parameters satisfy Equation 



( |12CI|) so that in front of the entire bracket there is an e-pole due to r(e) 



For every single / the double integral is finally calculated for e = numeri- 
cally and provides the coefficients of simple e-poles. Parts of the integrand 
with e 7^ lead to convergent contributions and are therefore omitted. 
The coefficient r(e)r(|) of the simple integral, however, contains an e^-pole. 
Hence, the integrand f{z = ^,w;e) must be expanded with respect to e up 
to the first order 



1 ^ r(f -e + a + 7 + /i-/3-2,5-2z/) 1 

2'"^''^ r(9-e-2z/-25 + /i + 7) ^ 4 



p. 1 ^p.l , r(g^ - 2{u -w)+ /x)r(^ - 2{6 + w)+^) 
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r(/3-i-e) ^ , ^ 2; 1 z 

— z— — ^ — ^ — V[b Vw^Viv 

r(4-e) ^ 2 ^ ^ 2 ^^2 2 ^^2 2 



2w - 7) + ^(2z/ - - - 2w - /x)) + 0(e^ 
2 



(125) 



where \l/(x) = denotes Euler's \1/ function and the relations ( |120|) and 
(|122|) have been used. The integral of the zeroth order in e provides the 
coefficients of the e^-poles and can even be executed analytically |^ 



27rz 4 ' ' 4 



w)T(- 



w)Ti- 



w] 



(126) 



T{5)T{5 + v~\)V{u)V{l 



T{6 



The integral of the first order in e yields the coefficients of the e-poles and 
is numerically calculated for every single I that belongs to this case 2 a), 
b) The whole integration is primitively convergent. 

This statement is true for rj_/(2, 4, 0, 1, 1, 2) and t_i_/(2, 4, 0, 1, 3, 3) in Equa- 
tion (|104|) . Here the parameters have the property 



2(/? + 5 + z/) - (a + 7 + /i) - (rf + 1) = 1 + e , 



(127) 



i.e. in front of the brackets in Equation (|124|) the coefficient is r(l + e) and 
consequently there is no e-pole. 

Therefore the double integral needs not to be computed, as it only leads 
to convergent contributions. The simple integral can be evaluated at e = 
because the coeffient contains only a simple e-pole. Due to the relation ( |122| ) 
that is also valid here, the simple integral is reduced to the one already solved 
in Equation ( |1 26| ) analytically. 

Now all integrals of type I from Equation (|104|) have been calculated. The 
way to solve the integrals of type F is step by step analogous to the method 
presented for type /. F(3; 1, 4, 0, 1, 1, 3) and F(l; 1, 3, 0, 1, 1, 3) from Equation 
( |104| ) belong to case 1, whereas F(l; 1, 4, 1, 2, 0, 1) belongs to case 2 a). 
Finally we present the sum of the two-loop diagrams for i up to convergent 
parts (cf. (pm) 



(1,1) 



1 



6^ lis 



0.01199e 



(128) 
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The 25 two-loop diagrams of the primitively divergent Fi 2 again only lead 
to integrals of type / and F 

X: ''^ = -4z^g|| [3/(1, 3, 0, 1, 1, 2) - /(I, 3, 1, 2, 2, 2) - 3/(3, 4, 0, 1, 1, 2) 
i=i P 

+ /(3, 4, 1, 2, 2, 2) - 2/(1, 3, 0, 1, 3, 3) + 2/(3, 4, 0, 1, 3, 3) 

- 3t^/(2, 4, 0, 1, 1, 2) + 2t^/(2, 4, 0, 1, 3, 3) + t±I{2, 4, 1, 2, 2, 2) 

- i/(0, 2, 1, 2, 1, 2) + /(2, 3, 1, 2, 1, 2) - i/(4, 4, 1, 2, 1, 2) 

- F(3; 1, 4, 0, 1, 1, 3) + F(3; 0, 3, 1, 2, 1, 3) + 2F(3; 1, 5, 1, 2, 0, 1)], 

(129) 

which can be evaluated with the methods demonstrated for Fi^i. This con- 
cludes the two-loop calculation. 
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Figure 1: Typical configurations of transverse (a) and longitudinal (b) or- 
dered phase. Regions of high density are shaded. 
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n+1 

(a) 





Figure 2: (a) Unsymmetric random potential, (b) random potential with 
equally high mountains, and (c) random potential with equally deep valleys. 
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Figure 4: The graphical elements of perturbation theory in the quasi-static 
model: (a) vertex, (b) Gaussian propagator, and (c) Gaussian correlator. 
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Figure 6: The complex integration path lying between two poles is decom- 
posed into two parts. 
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